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1. Introduction 

In SU (N c ) gauge theories without matter fields (gluodynamics) or matter in the adjoint representa- 
tion the action and measure are both invariant under center transformations while the trace of the 
Polyakov loop 



transforms non-trivially and therefore serves as an order parameter for the spontaneous breaking 
of center symmetric On the other hand the expectation value of the Polyakov loop is related to the 
free energie of an infinitely heavy quark, (P) oc exp(—f5 T F q ). At low tempertures gluons and quarks 
are confined in mesons and baryons and it needs an infinite amount of energy to free them, while at 
high temperatures confinement is lost and the energy becomes finite. Thus the Polyakov loop also 
serves as an order parameter for confinement and relates the breakdown of center symmetry to the 
confinement-deconfinement phase transition. In such gauge theories confinement is equivalent to 
the existence of an unbreakable string connecting a static quark anti-quark pair. 

In contrast in QCD or gauge theories with fundamental matter fields the center symmetry 
is explicitly broken and the Polyakov loop ceases to be an order parameter. As a consequence the 
string can break due to dynamical light quark production and in this sense confinement is equivalent 
to the existence of a string only at intermediate scales. It is widely believed that confinement is a 
property of pure gauge theories and to clarify the relevance of center symmetry it suggests itself to 
study pure gauge theories whose gauge groups have a trivial center. The exceptional Lie group 
is the smallest simple Lie group with this property which is simply connected. In a pioneering work 
the group in Bern has been convincingly demonstrated that G2 gluodynamics shows a first order 
finite temperature confinement-deconfinement phase transition [0, ||, §]. As in QCD confinement 
refers to confinement at intermediate scales, where a Casimir scaling of string tensions has already 
been reported [Q] and on large scales string breaking is expected to occur due to spontaneous gluon 
production []|] but so far not confirmed. 

Additionaly the gauge group SU (3) of strong interaction is a subgroup of G2 and this ob- 
servation has interesting consequences With a Higgs field in the fundamental 7 -dimensional 
representation one can break the G2 gauge symmetry to the SU (3) symmetry via the Higgs mech- 
anism. When the Higgs field in the action 



picks up a vacuum expectation value v, the 8 gluons belonging to SU (3) remain massless and the 
additional 6 gauge bosons acquire a mass proportional to v. In the limit v — > 00 they are removed 
from the spectrum such that G2 Yang-Mills-Higgs (YMH) theory reduces to SU (3) Yang-Mills 
theory. Even more interesting, for intermediate and large values of v the G2 YMH-theory mimics 
SU (3) gauge theory with dynamical 'scalar quarks'. 

The present paper deals with G2 gluodynamics in 3 dimensions and the G2 Gauge Higgs model 
in 4 dimensions. The simulations are performed with an efficient and fast implementation of a local 
hybrid Monte-Carlo algorithm. Below we shall calculate the potentials at intermediate scales for 
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Table 1: Representations of G2 with corresponding dimension and Casimir values. 

representation g [1,0] [0,1] [2,0] [1,1] [3,0] [0,2] [4,0] [2,1] 

dimension dgt 7 14 27 64 77 7T 182 189~ 

Casimir eigenvalue 12 24 28 42 48 60 72 64 

Casimir ratio 1 2 7/3 3.5 4 5 6 16/3 



static charges in the 8 lowest representations. We show that in 3 dimensions the string tensions 
on intermediate scales are proportional to the second order Casimir of the representations and for 
widely separated charges in the two fundamental representations we see a flattening of the potential 
which signals the breaking of the connecting string. In 4 dimensions we investigate the phase 
diagram a the G2 Gauge Higgs model and find a line of first order confinement-deconfinement 
phase transitions connecting G2 and SU (3) gluodynamics and two deconfmement phases separated 
by a second order phase transition. Details on the used algorithms and results can be found in [jg]. 

2. The group G 2 

G2 is the smallest of the five exceptional simple Lie groups and it is also the smallest simply con- 
nected simple Lie group with a trivial center. It has dimension 14 and rank 2. The fundamental 
representations are the defining 7 dimensional representation and the adjoint 14 dimensional repre- 
sentation. It is a subgroup of 50(7) and the gauge group SU (3) of strong interaction is a subgroup 
of G2. The corresponding coset space is a sphere [Q], 

G 2 /5£/(3)~5 6 , (2.1) 

meaning that every element % of G2 can be factorized as 

fy = ,9>.y with f€SU(3) and J^eG 2 /5£/(3). (2.2) 

This decomposition is used in our simulations to compute the exponential map of Q2 — > G2 [0]. 

Any irreducible representation of G2 is characterized by its highest weight vector /J. which 
is a linear combination of the fundamental weights, pL = p/J,m +?M(2)> w i m non-negative integer 
coefficients p,q called Dynkin labels. The dimension of an arbitrary irreducible representation 
S% = [p,q] can be calculated with the help of Weyl's dimension formula and is given by 

d S g = —(\+p)(\+q){2 + p + q){3+p + 2q)(4 + p + 3q){5 + 2p + 3q). (2.3) 

Below we also use the physics-convention and denote a representation by its dimension. For exam- 
ple, the fundamental representations are [1,0] =7 and [0, 1] = 14. An irreducible representation of 
G2 can also be characterized by the values of the two Casimir operators of degree 2 and 6. Below 
we shall need the values of the quadratic Casimir in a representation \p,q], given by 

«k = <€ M = 2p 2 + 6q 2 + 6pq + lOp + \8q. (2.4) 

For an easy comparison we normalize these 'raw' Casimir values with respect to the defining 
representation by Sf' = c to p . q /'io\.Q. The normalized Casimir values for the eight non-trivial repre- 
sentations with smallest dimensions are given in Tab. [I]. 
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Figure 1: Left panel: Polyakov loop expectation values at the finite temperature confinement-deconfinement 
transition on a 16 3 x 6 lattice. Right panel: Tunneling of the Polyakov loop and the (finite volume) phase 
coexistence at the first order transition at j3 cr it = 9.765. 



3. The confinement-deconfinement transition 

In SU (3) the situation is described as follows: Quarks and anti-quarks transform under the 
fundamental representations 3 and 3 and their charges can only be screened by particles with non- 
vanishing 3-ality, especially not by gluons. So in the confining phase the static quark anti-quark 
potential is linearly rising up to arbitrary long distances. As a consequence the free energy of a sin- 
gle quark gets infinite and the Polyakov loop expectation value vanishes. Hence in gluodynamics 
the Polyakov loop serves as order parameter for the Z3 centre symmetry and for confinement. 

In G2 we recall the decomposition of tensor products into irreducible representations, 

(7)®(7) = (1)©-", 

(3.1) 

(7)®(7)®(7) = (1)©--- 

The quarks in G2 transform under the 7-dimensional fundamental representation, gluons under the 
14-dimensional fundamental representation. Similarly as in SU (3) two or three quarks can build a 
colour singlet. Additionally three centre-blind dynamical gluons can screen the colour charge of a 
single quark, 

(7) ® (14) ® (14) ® (14) = (l)e-... (3.2) 

Thus the flux tube between two static quarks can break due to gluon production and the Polyakov 
loop does not vanish even in the confining phase. This shows that the Polyakov loop can at best 
be an approximate order parameter which changes rapidly at the phase transition and is small (but 
non-zero) in the confining phase. To characterise confinement we can no longer refer to a non- 
vanishing asymptotic string tension and vanishing Polyakov loop. Instead we define confinement 
as the absence of free colour charges in the physical spectrum. In the confining phase the static 
quark anti-quark potential rises linearly only at intermediate scales. Nevertheless we see a clear 
signal in the Polyakov loop and its distribution in the fundamental domain of G2 at the confinement- 
deconfinement transition (Fig. 
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4. Algorithmic considerations 
4.1 Local hybrid Monte-Carlo 



The corresponding lattice action for the G2 Yang-Mills-Higgs theory (1.2) reads 



5 Y mh[^,<J>] =j3£(l-itrRe^h) - K^x+fi%n^x, = 1, (4-1) 

where <J> is a seven component normalized real scalar field. Although there exists a heat-bath 
Monte-Carlo algorithm for G2 gluodynamics ^ we present a (local) HMC algorithm based on ^ 
for several good reasons: The formulation is given entirely in terms of Lie group and Lie algebra 
elements and there is no need to back-project onto G2, the autocorrelation time can be controlled 
(in certain ranges) by the integration time in the molecular dynamics part of the HMC algorithm 
and the inclusion of a (normalized) Higgs field is straightforward and does not suffer from a low 
Metropolis acceptance rate (even for large hopping parameters). The LHMC algorithm has been 
essential for obtaining the results in the present work. Since we developed the first implementation 
for G2 it is useful to explain the technical details for this exceptional group. As any (L)HMC 
algorithm for gauge theories it is based on a fictitious dynamics for the link- variables on the gauge 
group manifold. The "free evolution" on a semisimple group is the Riemannian geodesic motion 
with respect to the Cartan-Killing metric 

Without interaction the scalar field <J> is randomly distributed on the unit sphere and we set 

$ X = X $ O with G x eSO(l) (4.2) 
and constant 4>q. In the (L)HMC-dynamics the interaction term is given by the Yang-Mills-Higgs 



action ( |4.l[ ) of the underlying lattice gauge theory and in terms of the ^-variables we choose 



as Lagrangian for the HMC-dynamics 

L = -ij> {^x,^) 2 ~ ^£tr (^-') 2 - 5ymh[^, &l (4.3) 
^ x,n x 

where dot denotes the derivative with respect to the HMC time parameter t. The Lie algebra valued 
fictitious momenta conjugated to the link variable ^ )jU and site variable G x are given by 

m „ = — = -<% x u %- 1 , <p v = — — = -& x G-\ (4.4) 

The Legendre transform yields the following pseudo-Hamiltonian 

H = ~£tr*p^-i£tr^ + 5yMH[^,^]. (4-5) 

x,n x 

The equations of motion for the momenta are obtained by varying the Hamiltonian. The variation 
of Sym[% with respect to a fixed link variable yields the staple variable R x ^, the sum 
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of tiiple products of elementary link variables closing to a plaquette with the chosen link variable. 
Setting 

8y Xltl =%,ndT , 8%^ = ^dr = -^%^dT (4.6) 

with similar expressions for the momentum and field variables Sty x and 8ff x in the Higgs sector, 
one finds for the variation of the HMC-Hamiltonian 

8H = -Xtr^j^ -F^} -£tr^{$, -F x ) (4.7) 

X,fi X 

with the following forces in the gauge- and Higgs sectors 

F*,H = ^ [%,pRxji + K (^MI +Ii , Fx = ^v(L : v^ vPv) ' . (4.8) 

where the last sum extends over all nearest neighbors y of x and U xy denotes the parallel transporter 
from y to x. The variational principle implies that the projection of the terms between curly brackets 



in (4_/7) onto the Lie algebras 02 and so (7) vanish. Hence choosing a trace-orthonormal basis {T a } 
of 02 and {T/,} of so(7) the LHMC-equations read 

&x,n = - s $x,n%,n, &,M=I>(F XiM r fl )7; and 

(4.9) 

ff x = -% x ff x , % = ^ti(F x f h )f b . 

b 

The involved exponential maps are given by relatively simple analytical expressions ^ and a large 
step size (Leap frog second order integrator with 8x = 0.25, N t = 3 in most of our simulations) 
allows for a fast and efficient implementation of the algorithm. 

4.2 Exponential error reduction for Wilson loops 

In the confining phase the rectangular Wilson loop scales as W(L,T) oc exp(— oh - T). In order 
to estimate the string tension a we probe areas LT ranging from up to 100 and thus W will 
vary by approximately 40 orders of magnitude. A brute force approach where statistical errors 
for the expectation value of Wilson or Polyakov loops decrease with the inverse square root of 
the number of statistically independent configurations by just increasing the number of generated 
configurations will miserably fail. Thus to obtain accurate and reliable numbers for the static 
potential and to detect string breaking we implemented the multi-step Liischer-Weisz algorithm 
with exponential error reduction for the time transporters of the Wilson loops With this method 
the absolute errors of Wilson lines decrease exponentially with the temporal extent T of the line. 
This is achieved by subdividing the lattice into n t sublattices V\,... ,V nt containing the Wilson loop 
and separated by time slices plus the remaining sublattice, denoted by V, see Fig. |2| (left panel). 
At the first level in a two-level algorithm the time extent of each sublattice V n is 4 such that n t is 
the smallest natural number with 4n t >T + 2. In Fig. || (left panel) T = 14 and the lattice is split 
into four sublattices Vi,V2,V3,V4 containing the Wilson loop plus the complement V. The Wilson 
loop is the product of parallel transporters W = ^'^'^ £^,^2^1- If a sublattice V„ contains only 
one connected piece of the Wilson loop (as V\ and V4 do) then one needs to calculate the sublattice 
expectation value 

(%)n = ^r I ®<%^ n e- s , (4.10) 

sublattice n 
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Figure 2: (Color online) Iterative slicing (from left to right) of lattice and Wilson loop during the multi-level 
algorithm. 

if V n contains two connected pieces (as V2 and V3) then one needs to calculate <X> ^') H . The 
updates in each sublattice are done with fixed link variables on the time-slices bounding the sublat- 
tice. Calculating the expectation value of the full Wilson loop reduces to averaging over the links 
in the n t + 1 time slices, 

(w) = (^((^ 1 ) 1 (^ 2 »^ 2 / ) 2 ---(^ t _ 1 »^;_ 1 ), Jt _ 1 (5; t ), lt )) . (4.ii) 

\ \ / / boundaries 

Here ^ is that particular contraction of indices that leads to the trace of the Wilson loop. In a 
two-level algorithm each sublattice V n is further divided into two sublattices V n \ and V, 12 , see 
Fig. U (middle panel), and the sublattice updates are done on the small sublattices V n> k with fixed 
link variables on the time slices separating the sublattices V n 1. This way one finds two levels of 
nested averages. Iterating this procedure gives the multi-level algorithm. Since the dimensions dgg 
grow rapidly with the Dynkin labels [p,q] - for example, below we shall verify Casimir scaling for 
charges in the 189 dimensional representation [2,1] - it is difficult to store the many expectation 
values of tensor products of parallel transporters. Thus we implemented a slight modification of the 
Liischer-Weisz algorithm where the lattice is further split by spatially slicing along a hyperplane 
orthogonal to the plane defined by the Wilson loop, see Fig. § (right panel). 

In the present work we use a two-level algorithm with time slices of length 4 on the first 
and length 2 on the second level to calculate (W) for Wilson loops (and hence transporters J^ n ) 
of varying sizes and in different representations and a three-level algorithm with times lices 8,4 
and 2 for Polyakov loops. To avoid the storage of tensor products of large representations we 
implemented the modified algorithm as explained above. 



5. String tension and Casimir scaling in G% gluodynamics 

The static inter-quark potential is linearly rising on intermediate distances and the corresponding 
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string tension will depend on the representation of the static charges. We expect to find Casimir 
scaling where the string tensions for different representations M and St' scale according to 

°— (5.1) 

cm c%> 

with quadratic Casimir c@. Although all string tensions will vanish at asymptotic scales it is still 
possible to check for Casimir scaling at intermediate scales where the linearity of the inter-quark 
potential is nearly fulfilled. Up to these scales we can parametrize the potential with 

V*(R) = Y*-^ + osR. (5.2) 

K 

To extract the static quark anti-quark potential two different methods are available. The first 
makes use of the behavior of rectangular Wilson loops in representation ffl for large T, 

(WM(R,T))=exp(K M (R)-V M (R)T). (5.3) 

The potential can be extracted from the ratio of two Wilson loops with different time extent accord- 
ing to 

' T (Ws,(R,T + x)) 

We calculated the expectation values of Wilson loops with the two-level Luscher-Weisz algorithm 



and fitted the right hand side of ( |5.4] ) with the potential Vm(R) in (p^2|). The fitting has been done 
for external charges separated by one lattice unit up to separations 7? with acceptable signal to noise 
ratios. From the fits we extracted the constants y^, CCgg and a@ entering the static potential. For an 
easier comparison of the numerical results on lattices of different size and for different values of /3 
we subtracted the constant contribution to the potentials and plotted 

V®(R)=V m {R)-y M (5.5) 

in the figures. The statistical errors are determined with the Jackknife method. In addition we 
determined the local string tension 



<7loc,^ [R + 2 ) = ' ( 5 - 6 ) 



given by the Creutz ratio 



The second method to calculate the string tensions uses correlators of two Polyakov loops, 

V 3g {R) = -^-\n{P^)P^{R)). (5.8) 
Pt 

The correlators are calculated with the three-level Luscher-Weisz algorithm and are fitted with the 
static potential V^(R) with fit parameters y@,CC& and o^g. Now the local string tension takes the 
form 

p\ 1 {Pm(0)Pm(R+P)) 



2 J p T p (Pa{0)P<g{R)) 
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Figure 3: (Color online) Left panel: Continuum scaling of the fundamental potential. Right panel: Unsealed 
potential with /3 = 40 on a 28 3 lattice for different representations. 



-1.0 



11 = 7 
11 = 14 
11 = 27 
11 = 64 

n = n 

K = 77' 
7? = 182 

re = 189 



0.0 



0.5 



1.0 1.5 



2.0 



a T (B) 



2.5 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
fiR 



Figure 4: (Color online) Left panel: Scaled potential with /3 = 40 on a 28 3 lattice. Right panel: Ratio of 
the local string tension with j3 = 40 scaled on a 28 3 lattice for the eight smallest representations. 



5.1 Casimir scaling in 3 dimensions 

Most LHMC simulations are performed on a 28 3 lattice with Wilson loops of time extent T = 12. 



To extract the static potentials from the ratio of Wilson loops in (5.4) we chose z = 2. To check for 
scaling we plotted the potentials in 'physical' units, V//i, with mass scale set by the string tension 
in the 7 representation, 

M = ^, (5-10) 

as function of /J.R in Fig. || (Left panel). We observe that the potentials for the three values of j8 are 
the same within error bars. In addition they agree with the potential (in physical units) extracted 
from the Polyakov loop on a much larger 48 3 lattice. 

In Fig. H (Right panel) we plotted the values for the eight potentials Vj, ■ . ■ , Vi89 (with statistical 
errors) measured in 'physical units' jj. defined in ( 5.10| ). The distance of the charges is measured in 



the same system of units. The linear rise at intermediate scales is clearly visible, even for charges 
in the 189 dimensional representation. 
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Fig. |] (Left panel) contains the same data points rescaled with the quadratic Casimirs of the 
corresponding representations. The eight rescaled potentials fall on top of each other within error 
bars. This implies that the full potentials for short and intermediate separations of the static charges 
show Casimir scaling. 

To further check for Casimir scaling we calculated the local string tensions with p = 1 for all 
R between 1 and 10. The horizontal lines in Fig. |] (Right panel) are the values predicted by the 
Casimir scaling hypothesis. Clearly we see no sign of Casimir scaling violation on a 28 3 lattice 
near the continuum at /3 = 40. Of course, for widely separated charges in higher dimensional repre- 
sentations the error bars are not negligible even for an algorithm with exponential error reduction. 



5.2 String breaking and glue-lumps in 3 dimensions 

To observe the breaking of strings connecting static charges at intermediate scales when one further 
increases the separation of the charges we performed high statistics LHMC simulations on a 48 3 
lattice with j8 = 30. We calculated expectation values of Wilson loops and products of Polyakov 
loops for charges in the two fundamental representations of Gi. When a string breaks then each 
static charge in the representation M at the end of the string is screened by N(M) gluons to form 
a color blind glue-lump. We expect that the dominant decay channel for an over-stretched string 
is string — > gluelump + gluelump. For a string to decay the energy stored in the string must be 
sufficient to produce two glue-lumps. According to ( |3.2| ) it requires at least 3 gluons to screen a 
static charge in the 7 representation, one gluon to screen a charge in the 14 representation and two 
gluons to screen a charge in the 27 representation. We shall calculate the separations of the charges 
where string breaking sets in and the masses of the produced glue-lumps. The mass of such a 
quark-gluon bound state can be obtained from the correlation function according to 



exp(-m,^T)ocC^(T) 



F nv(y) 



n=\ 



(5.11) 



:»'.b i 



where M{% yx ) is the temporal parallel transporter in the representation ffl from x to y of length T. 
It represents the static sources in the representation M. The vertical line means projection of the 
tensor product onto that linear subspace on which the irreducible representation M acts, 



(14® 14( 



14) =StK 



(5.12) 



For example, for charges in the 14 representation the projection is simply 



14,a 



where F£ v T a 



(5.13) 



For charges in the 7 representation we must project the reducible representation 14(g) 14(g) 14 onto 
the irreducible representation 7. Using the embedding of G2 into 50(7) representations one shows 
that this projection can be done with the help of the totally antisymmetric £-tensor with 7 indices, 



Ffiv (x) ® Fp V (x) (g) F MV (x) 



T.a 



FP v (x)F^(x)F^ v (x)e ahcdef J b % q e T^. (5.14) 



Fig. U (Left panel) shows the logarithm of the glue-lump correlator (5.11) as function of the sepa- 
ration of the two lumps for static charges in the fundamental representations 7 and 14. The linear 
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Figure 5: (Color online) Left panel: Glue-lump correlator (lattice size 48 3 , j3 = 30). Right panel: Potential 
for both fundamental representations at j3 = 30 and corresponding glue-lump mass 
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Figure 6: (Color online) Left panel: Local string tension (48 3 lattice, /3 = 30). Right panel: Casimir scaling 
of local string tension (48 3 lattice, j3 = 30). 



fits to the data yield the glue-lump masses 

m 7 a = 0.46(4), m l4 a = 0.761(3). (5.15) 
Thus we expect that the subtracted static potentials approach the asymptotic values 

%g — >2m<%-y<%. (5.16) 
With the fit-values 77a = 0.197(1) and y u a = 0.381(2) we find 

V 7 /jU — ^ 3.47 , Vu/n — ^5.47. (5.17) 

Fig. H (Right panel) shows the rescaled potentials for charges in the fundamental representations 
together with the asymptotic values (5.17) extracted from the glue-lump correlators. At fixed cou- 
pling jS = 30 both potentials flatten exactly at separations of the charges where the energy stored 
in the flux tube is twice the glue-lump energy. 
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Fig. |6] shows the local string tensions in the two fundamental representations and their ratio. 
Although more gluons (three instead of one) are involved we see clearly that the string connecting 
charges in the adjoint representation breaks earlier than the string connecting charges in the defining 
representation. 



6. 4-dimensional Gauge-Higgs model 



The lattice action (41) depends on the inverse gauge coupling j3 and the hopping parameter 
K which is proportional to the vacuum expectation value of the Higgs field. For j3 — > °° the gauge 
bosons decouple and the theory reduces to an SO (7) -invariant nonlinear a-model which shows 
spontaneous symmetry breaking down to SO(6) at some critical value fQ.. The phase transition is 
of second order. 



For k — > oo the factor ,y in the decomposition (2^2) is frozen and we end up with an SU (3) 
gauge theory for the factor which shows a weak first order deconfinment transition. With respect 
to the unbroken subgroup SU(3) the fundamental representations (7) and (14) branch into the 
following irreducible SU (3) -representations: 

(7)— ► (3) ©(3)0(1) 
(14)^(8)0(3)0(3). (6.1) 

The Higgs field branches into a scalar quark, scalar anti-quark and singlet with respect to SU (3). 
Similarly, a G2-gluon branches into a massless SU (3)-gluon and additional gauge bosons with 
respect to SU(3). The latter eat up the the non-singlet scalar fields such that the spectrum in 
the broken phase consists of 8 massless gluons, 6 massive gauge bosons and one massive Higgs 
particle. If K is lowered, in addition to the 8 gluons of SU(3), the 6 additional gauge bosons of 
G2 with decreasing mass begin to participate in the dynamics. Similarly as dynamical quarks and 
anit-quarks in QCD, they transform in the representations (3) and (3) of SU (3) and thus explicitly 
break the Z3 center symmetry. As in QCD they are expected to weaken the deconfinement phase 
transition. Thus it has been conjectured in [Q] that there may exist a critical endpoint where the 
deconfinement transition disappears. For K = we recover G2 gluodynamics with a first order 
deconfinement phase transition. 

We measure the Polyakov loop as an (approximate) order parameter for confinement and in- 
vestigate the corresponding critical curve in the /3-fc plane. For large K the confinement phase 
in St/ (3) is characterised by (xi) = 1. If in the deconfinement phase the Z 3 centre symmetry 
of the remaining SU (3) is broken then the ambiguity of measuring * s fixed by choosing the 
Polyakov loop that points into the 1 -direction of St/ (3). Technically this is achieved by taking 
(Xi) — > 3 — 2(xi)- The results are shown in Fig. ^ (Left panel). We observe that on the small 
lattice the Polyakov loop jumps along a continuous line connecting the deconfinement transitions 
of pure G2 and pure SU (3) gluodynamics. This points to a continuous line of deconfinement tran- 
sitions all the way from fc = to fc = 00. To see whether this is indeed the case we performed 
high-precision simulations on larger lattices. A careful analyses of histograms and suszeptibilies 
for Polyakov loops and the Higgs part of the action confirm the results on the small lattice. Unfor- 
tunately a rather small region in parameter space is left where we cannot resolve the order of the 
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Figure 7: Left panel: Phase diagram of the Gauge-Higgs model in terms of the Polyakov loop expectation 
value (xi) on a 12 3 x 2 lattice. Right panel: Phase diagram of the Gauge-Higgs model on a lattice with 6 
links in temporal direction (mainly computed on a 16 3 x 6 lattice). Note the different scales in the separate 
regions of the diagram. First and second order (or crossover) transitions are marked with a continuous and 
dashed line (to guide the eye), respectively. 



transition. In this small region there could exist a crossover from the confining to the deconfining 
phase. 

On a larger (16 3 x 6) lattice we calculate also the full phase diagram including the Higgs 
S0{1) — > SO(6) transition (Fig. [7], Right panel). Phase transitions are obtained by observing 
susceptibility peaks in the Polyakov loop and the Higgs part of the action. We also investi- 
gate the order of the confinement-deconfinement transition using the histogram method for the 
Polyakov loops. For the order of the Higgs transition line we consider the finite size scaling of 
3k(V X "Lx,^x+jiUx,ii^x) (n = 1,2) for lattices up to 20 3 x 6. With the data obtained so far the 
point where the second order SO(J) — > SO(6) transition may turn into a crossover cannot be de- 
termined reliably. If the triple point exists then an extrapolation to the point where the confining 
phase meets both deconfining phases leads to the couplings j8 tr ip = 9.62(1) and fc tr j p = 1.455(5). 



7. Conclusions 

In the present work we implemented an efficient and fast LHMC algorithm to simulate G2 gauge 
theory in three and four dimensions. In addition we implemented a slightly modified Liischer- 
Weisz multi-step algorithm with exponential error reduction to measure the static potentials for 
charges in various G2 representations. The accurate results in 3 dimensions show that the static 
potentials show Casimir scaling on intermediate scales within a few percent statistical errors. Thus 
we conclude that in 3 dimensional G2 gluodynamics the Casimir scaling violations of the string 
tensions are small for all charges in the representations with dimensions 7, 14,27,64,77,77', 182 
and 189. 

For larger separations we detect string breaking in the two fundamental representations exactly 
at the expected scale where the energy stored in the flux tube is sufficient to create two glue-lumps. 
To confirm this expectation we calculated masses of glue-lumps associated with static charges in 
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the fundamental representations. Here, close to the string breaking distance, systematic Casimir 
scaling violations show up (for a more detailed discussion see 

In 4 dimensions we explored the full phase diagram of the G2 Gauge Higgs model. We found a 
line of first order phase transitions connecting G2 and SU (3) gluodynamics and a line of second or- 
der phase transitions separating the G2 deconfinement phase from the SU (3) deconfinement phase. 
Unfortunately we cannot exclude the existance of a small crossover region near the would-be tripel 
point of the system. Details on the G2 Gauge Higgs model will be published in a follow up paper. 
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